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ABSTRACT 


Linearized  free -streamline  theory  is  used  to  calculate  the  spray- 
sheet  thickness  and  lift-slope  for  a  flat  plate,  cavitating,  two-dimensional 
hydrofoil  in  a  channel  of  finite  depth  with  an  upper  free  surface  and  lower 
boundary  partly  free  and  partly  rigid.  Only  the  case  of  zero  cavitation 
number  is  considered.  Seme  measurements  were  made  of  the  submergence 
of  a  hydrofoil  of  four  inch.; 3  chord  beneath  the  undisturbed  free  surface  at 
velocities  of  12  and  18  ft.  /sec.  These  agree  with  the  trends  of  the  theory 
but  not  with  the  magnitudes,  the  submergence  always  being  greater  than 
that  predicted  by  the  theory. 
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1.  Introduction 


Recently,  several  papers  dealing  with  the  presenc  e  of  a  free 
surface  in  the  proximity  of  a  lifting,  fully  cavitating  hydrofoil  have  appear¬ 
ed.  Of  these,  one  by  Dawson  and  Bate^  ^  is  primarily  experimental. 

In  their  work  a  cavitating  or  ventilated  flat  plate  hydrofoil  was  investigated 
near  a  free  surface.  The  submergence  was  varied  so  that  the  hydrofoil 
went  from  a  cavitating  condition  to  a  planing  condition  near  the  surface. 
Forces  were  measured  in  this  experiment  and  compared  to  the  results  of 
several  free-streamline  theories  applicable  to  this  type  of  flow.  One  of 
these  is  a  linearized  theory  for  the  flow  of  an  infinitely  deep  stream  past 
a  cavitating  hydrofoil  in  the  presence  of  a  free  surface.  This  is  the  linear¬ 
ized  theory  due  to  Schot.  ^  Johnson  ^  considers  both  flat  plate  and  cam¬ 
bered  hydrofoils  of  finite  aspect  ratio  near  a  free  surface.  This  latter 
work  is  largely  experimental.  In  these  cases  mentioned,  both  experiment 
and  theory  are  concerned  with  hydrofoils  at  small  angles  of  attack  relative¬ 
ly  near  a  free  surface,  i.  e.  two  chords  and  less,  and  with  small  cavitation 
numbers  approaching  zero.  Furthermore,  the  influence  of  gravity  has  been 

everywhere  neglected.  The  most  complete  vaeory  which  has  so  far  been 

(4) 

produced,  is  ar.  exact  non-linear  theory  due  to  Green.  He  considers  in 
one  of  his  papers  the  planing  of  a  two-dimensional  flat  plate,  without  the 
influence  of  gravity,  on  a  stream  of  finite  depth.  Although  Green's  work  is 
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exact  and  the  linear  theories  reported  by  Schot  and  referred  to  by  Daws  or 
and  Bate  are  not#  Green,  however,  did  not  compute  a  sufficient  number 
of  cases  of  interest  to  users  of  water  tunnels.  Furthermore,  his  solutions 
are  not  presented  in  an  easily  usable  form.  All  free  surface  water  tunnel 
experiments  are  conducted  in  a  channel  of  finite  depth,  as  opposed  to  the 
analyses  of  Refs.  (2)  and  (3).  The  measurements  of  Dawson  and  Bate, 
particularly  near  the  free  surface,  do  not  agree  well  with  the  calculations 
of  the  submergence  effect  to  be  found,  for  example,  in  Ref.  (2),  and  it  was 
suggested  in  Ref.  (1)  that  the  discrepancies  are  probably  not  due  only  to 
the  influence  of  gravity,  at  least  as  estimated  by  current  planing  theories. '  ' 
Briefly,  the  submergence  effect  on  the  lift  of  a  flat  plate  hydrofoil  is  to  in¬ 
crease  the  lift  force  as  the  Surface  is  approached.  In  the  limiting  case  of 
zero  cavitation  number,  the  cavity  is,  of  course,  infinitely  long  and  the 
lift-slope  goes  from  the  value  of  x/2  for  infinitely  deep  submergence  to 
the  value  ir  at  the  surface  for  infinitesimal  angles  of  attack-  In  the  experi¬ 
ment  of  Ref.  (1)  the  lift  increases  more  slowly  as  the  surface  is  approached 
than  is  indicated  by  linearized  theory  for  flows  of  infinite  depth  and  zero 
cavitation  number.  Incidentally,  the  lift-slopes  are  reasonably  well  ac¬ 
counted  for  by  the  theory.  Until  a  complete  solution  accounting  for  all  of  the 
effects  of  non-zero  cavitation  number,  channel  depth  and  neighboring  free 
surface  is  available,  the  exact  cause  of  this  discrepancy  probably  cannot 
be  determined.  In  the  present  report,  however,  it  has  been  decided  to 
remove  at  least  one  of  the  drawbacks  of  the  theoretical  work  mentioned 
previously;  namely,  there  is  a  bottom  at  a  finite  distance  beneath  the  hydro- 
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foil.  Of  course,  Green's  exact  theory  accounts  for  this  effect.  It  was  be¬ 
lieved,  however,  that  it  was  no  more  work  to  solve  the  problem  anew  by 
means  of  the  linearized  free  streamline  theory  than  to  recalculate  Green's 
exact  solution.  In  fact,  because  of  the  flexibility  afforded  by  the  linear 
theory  in  determining  the  flow  past  relatively  complex  shapes  and  non¬ 
zero  cavitation  numbers,  it  v/as  decided  to  follow  this  course  rather  than 
that  of  the  exact  theory.  The  linear  theory,  while  suffering  notable  limita¬ 
tions  is  expected  to  provide  a  good  value  for  the  lift-slope  at  vanishing 
angles  of  attack  and  to  account  properly  for  the  effects  of  submergence 
depth  of  channel,  and  even  for  the  effect  of  a  finite  length  of  the  bottom 
channel  if  such  is  the  case. 

The  results  of  the  calculations  to  be  presented  herein  should  be 
of  some  value  to  the  users  of  water  tunnels.  Extensive  tables  and  charts 
of  the  various  force  and  geometric  quantities  involved  are  presented.  In 
addition,  the  method  of  computing  the  results  and  the  approximations  for 
digital  computations  used  herewith,  are  discussed  in  some  detail  so  that 
the  reader  may  calculate  additional  cases  for  himself.  In  a  subsequent 
report  it  is  planned  to  consider  the  case  of  non-zero  cavitation  numbers. 

2.  Formulation  of  the  problem 

We  shall  apply  the  free  streamline  theory  to  the  steady,  two-dim¬ 
ensional  flow  past  a  flat  plate.  The  stream  is  of  finite  depth  with  the  upper 
surface  bounded  by  a  free  surface,  and  the  lower  surface  bounded  partly  by 

A  companion  report  (Ref.  (8))  may  also  be  consulted  for  the  computaion- 
al  method  used. 
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a  solid  wall  and  partly  by  a  free  jet.  The  plate,  which  has  both  sharp 
leading  and  trailing  edges,  is  inclined  at  a  certain  angle  of  incidence.  A 
spray  sheet  is  formed  at  the  leading  edge  but  the  flow  off  the  trailing  edge 
is  smooth.  The  sharpness  of  the  edges  allows  one  to  expect  that  the 
separations  of  flow  occur  at  these  points  as  shown  in  Fig.  1.  Far  down¬ 
stream,  the  spray  sheet  and  the  free  jet  become  inclined  at  certain  angles 
to  the  oncoming  flow.  The  magnitude  of  the  velocity  along  these  surfaces 
equals  that  of  the  free  stream. 

Certain  simplifications  are  then  introduced  into  the  problem: 

i)  the  gravity  effect  is  neglected  throughout  the  analysis, 

ii)  the  cavitation  number  is  zero  (infinitely  long  cavity), 

iii)  the  angle  of  incidence  of  the  plate  with  the  free  stream 
as  well  as  the  final  inclinations  of  the  spray  sheet  and 
jet  are  assumed  to  be  small  so  that  linearization  of  the 
problem  is  permitted. 

If  we  denote  the  velocity  by  q  then  a  complex  perturbation  velo¬ 
city  function 


w(z)  =  u  -  iv  (2.1) 

can  be  introduced  such  that  q  =  U(1  +  u,  v). 

The  dimensionless  perturbation  velocity  components  u(x,y)  and 
v(x,y)  are  assumed  to  be  much  smaller  than  unity.  Applying  the  usual 
linearized  boundary  condition  on  the  free  surfaces,  the  plate  and  the  solid 
wall,  we  have 
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u  =  0  on  all  free  surfaces  (constant  pressure  and  speed) 
v  =  -a  on  the  wetted  side  of  the  flat  plate, 
v  =  0  on  the  solid  wall  forming  part  of  the  lower  boundary. 

The  linearized  version  of  the  above  problem  may  then  be  depicted 
schematically  in  a  complex  z-plane  as  shown  in  Fig.  2. 

One  would  notice  that  after  linearization,  the  stagnation  point  on 
the  plate  coincides  with  the  leading  edge,  furthermore,  the  streamline 
at  that  point  splits  into  two  branches. 

In  addition  to  the  above  boundary  conditions,  the  perturbation 
velocity  components  must  satisfy  the  following  physical  conditions: 

i)  there  is  a  singularity  at  the  leading  edge, 

ii)  the  streamline  has  a  continuous  slope  at  the  plate- 
cavity  juction  (Kutta  condition), 

iii)  the  perturbation  velocity  components  vanish  upstream 
(at  x  =  -  oo ). 

3.  Preliminary  results 

We  denote  the  depth  of  the  stream  at  upstream  infinity  by  H,  the 
thickness  of  the  spray  sheet  by  6,  and  the  final  inclinations  of  the  spray 
and  jet  by  6  and  6  respectively.  These  quantities ,  essential  in  our 
calculations  of  the  force  coefficient  and  the  flow  geometry,  are  to  be  deter¬ 
mined  from  the  theory.  In  addition,  the  following  quantities  are  specified: 
U  the  velocity  of  the  free  stream,  the  length  of  the  plate,  a  the 

angle  of  incidence,  h  the  height  of  the  plate  leading  edge  above  the  bottom 
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wall,  and  I  the  length  of  the  bottom  solid  wall  downstream  of  the  leading 
edge.  The  submergence  of  the  hydrofoil  beneath  the  undisturbed  oncoming 
flow  is  the  difference  between  H  and  h,  i.  e.  ,  H  -  h.  We  shall  now  make 
a  preliminary  calculation  based  on  continuity  and  momentum  considerations 
which  gives  some  useful  relations  between  the  forces  and  flow  geometry. 

The  ultimate  thickness  of  the  spray  sheet  and  lower  jet  must 
equal  the  upstream  depth  by  the  continuity  relation.  Thus  if  t  =6/ H,  the 
thickness  of  the  jet  is  (1  -  «  )H.  The  drag  on  the  plate  can  be  found  from 
a  momentum  balance  and  in  terms  of  c  it  is 

D  =  ^  ^  6*  +  (1  (3.1) 

using  the  approximation  cos©  =  1  -  0l/2  etc.  When  the  hydrofoil  is  a 
flat  plate  the  force  on  the  plate  must  be  normal  to  the  surface.  Thus,  if 
L  and  D  are  the  lift  and  drag  force  perpendicular  to  and  parallel  to  res¬ 
pectively  the  oncoming  flow  of  speed  U, 

D  =  L  tan  a  =  L  a  (3.2) 

where  a,  assumed  small,  is  the  angle  of  attack  of  the  plate  with  respect 
to  the  oncoming  flow. 

The  lift  coefficient  is  defined  as 


hence 
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Given  the  depth  above  the  bottom,  h,  and  the  upstream  depth, 
H,  we  may  expect  that  the  spray  sheet  thickness,  6,  and  hence  <  will 
vary  with  angle  of  attack.  We  anticipate  that  in  the  linear  theory  already 
sketched  in  section  (2),  velocities  and  hence  angles  such  as  6 1  ,  etc.  are 
proportional  to  a.  the  angle  of  attack.  Thus 


C,  = 


6a 


and  the  ratios  0^/a  etc.  ,  will  be  independent  of  a  in  the  present  theory. 
But  6,  within  our  linear  theory,  can  contain  a  linear  term  in  a.  Now 
in  the  limit  as  a  approaches  zero,  the  spray  sheet  thickness  w'ill  ap¬ 
proach  d  .  Thus  we  can  put 

6  =  d  ~a  g  (3.4) 

where  g  is  a  function  only  of  the  geometrical  ratios  C^/d,  C^/h  which 
has  to  be  determined  from  the  theory.  For  vanishingly  small  angles  of 
attack  we  will  have 


and  for  larger  angles  (restricting  ourselves  now  to  the  case  of  a  rigid, 

bottom  for  which  0  is  zero  )  we  have 

2 


(3.  5) 

flat 


(3.6) 


We  will  now  develop  the  linear  theory  to  determine  the  quantities 

(d/a),  (6  /a),  and  g.  With  these  and  given  values  of  d.  h.  and 

1  2 

the  spray  sheet  thickness,  submergence  and  lift  of  the  hydrofoil  can  be 


determined. 
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4.  Theory 

Consider  the  transformation 


=  -  b  log  (1  +-^)  -  c  log  (1 


(4.  1) 


where  £  -  4  +  iq. 

The  flow  in  the  physical  plane  is  mapped  into  the  upper  half  t,  -plane 

and  the  entire  boundary  of  the  flow  goes  into  the  real  axis,  shown  in  Fig.  3. 

The  leading  and  trailing  edges  of  the  plate  in  the  z-plane  are  at  the  origin 

and  -  1  respectively.  Upstream  infinity,  A,  transforms  into  infinity.  The 

« 

end  of  the  solid  wall,  B  ,  the  free  jet,  B,  and  the  spray  sheet,  C,  have  their 

I 

corresponding  image  points  at  -  b  ,  b  and  c.  The  value  of  b  approaches 

•  i 

b  as  the  length  I  is  extended  to  infinity.  The  quantities  b  ,  b  and  c  are 
parameters  whose  numerical  values  will  be  determined  later  from  the  map- 
ping  function.  We  know,  however,  at  this  stage  c  >  0,  b  >  b>  1. 

The  complex  perturbation  velocity  in  £-plane,  denoted  by 


v  ( £  )  =  u  -  iv 


takes  the  following  boundary  values  on  the  Rf  4  or  4  axis. 

v  =  0  -  oo  <  ^  <  -  b 


t 


u  =  0 

-  b  <  4  < 

- 1 

v  =  -a 

-  *  <  4  < 

0 

u  =  0 

o  <  4  < 

00 

The  flow  in  the  upper  half  £  -plane  can  be  continued  analytically 
into  the  lower  half  by 


v(£)  -  -v(&)  . 


(4  2) 


On  the  Rl  £  axis,  two  branch  cuts  appear;  one  starts  from  £  -  -  oo 

I 

to  -b  and  the  other  from  |  =  -  1  to  the  origin.  These  cuts  correspond 
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to  the  plate  and  the  solid  wall  in  the  physical  plane  respectively.  We  now 


> 


ne 


■s 


re 


have  a  boundary  value  problem  involving  the  unknown  function  v  (£ )  which 
can  be  formulated  and  solved  as  a  Hilbert  problem. ^ 

We  introduce  the  auxiliary  function 


H(U  -  i 


which  as  the  following  properties 


vtnrRT-r 


bT’ 


(4.  3) 


i)  H(4 )  has  the  proper  branch  cuts  (t,  -  0,  -  1 ,  -  b'  and  oo 
are  branch  points). 

ii)  It  provides  the  correct  leading  edge  singularity  (square 
root  singulaity). 


iii)  On  the  Ri  t,  axis; 
Im  H  =  0 
Rf  H  =  0 
Rf 

Ri  H  -  0 


H  -  / 

H  '  VTT7 


4  )(b’  +  4 ) 


-oo<  4  <  -b' 
_b’<  4  <  -l 
Im  H  =  0  -1  <  4  <  0 

0  <  4 


A  new  function 


G(4)  =  v(C,)  H(4) 
is  then  formed  which  has  on  the  Rf  4  axis; 
ImG  =  0 
Im  G  ~  0 
Im  G  =  -a 
ImG  -  0 


(4.  4) 


'-4 

0  +  4)(b 1 +  4T 


— oc  <4  <  —  b' 

-b’<  4  <  -i 
-l  <  4  <  o 
o  <  4  • 


(6) 


With  the  aid  of  Plemelj’s  formula.  the  solution  of  the  Hilbert 
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problem  is  given  immediately  by 


& 


(4.5) 


In  order  to  evaluate  the  final  inclinations  0  and  0  we  notice 

l  2 

that  according  to  the  linearized  theory,  the  streamline  has  a  slope 


a£  =  e  = 


(4.6) 


In  addition,  the  boundary  conditions  on  the  free  surface  require 


that  everywhere  u  =  0.  Hence,  we  obtain 

o 


and 


®»  =  -^v(-b>=  :bl-J  pti  \fuT 0(1'+$  dt ' 


(4.7) 

(4.8) 


The  mapping  function  gives  us  the  following  relations  of  the  de¬ 
sired  parameters. 

i  l 

(4.  9) 


ch=  -  b  log(l  "  ~  c  logU  +  “  )  * 


B'  =  -b  logO  +  £■  )-c  log(l  -  ^)  , 


(4.  10) 


where  B'  =  1  -  i it  h  , 

and  h  =  it  b  .  (4.11) 

By  specifying  the  ratio  C^/h  and  the  length  of  the  wall  I  ,  values  of 
b1,  b  and  c  can  be  calculated  from  these  simultaneous  equations.  The 
detailed  evaluation  of  these  quantities  will  be  given  in  a  later  part  of  the 
report. 

The  geometry  of  the  problem  can  also  be  solved  since  the  velocity 


No  complimentary  function  is  needed  here  and  Eq.  (4.  5)  is  the  complete 
solution. 
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field  is  now  known.  To  simplify  the  matter,  the  limiting  case  of  an  infinite¬ 
ly  long  solid  wall  will  be  worked  out  (  for  which  b  b').  The  reader  may 

refer  all  quantities  in  the  following  calculation  to  Fig.  4. 

o 

( 

Denote 


thus 


J  =  J,t  -  c  \Z(l  +  t)(b+t)  • 

e- 

1  *  y  < 


and 


)(b+c ) 


6  =  0 
2 

_  6  {c+1  )(c+b) 


(4.  12) 

(4.  13) 
(4.14) 
(4.  15) 


Since  b>  1  and  c  >  0  at  all  times  we  may  introduce  new  para¬ 
meters  k  and  k'  such  that 

k2  =  1 

D 

and  k'2  =  1  -  k2  <  1 

(7 ) 

The  integral  J  can  now  be  written  in  terms  of  known  eliptic  functions 
in  the  form 

J  =  2kijl  ,k)  -  K]  (4.  16a) 


or  in  an  alternate  form 


J  =  1I\/(UcHbtc)  Ao  -  Irf  k 


(4.  16b) 


where 


AQ(i|i,k)  =  ~[  EF(^.k’)  +  KE(4>,  k* )  -  KF(4>»  k' )] 


is  Heuman's  Lambda  function,  K  and  E  are  the  complete  elliptic  integrals 
of  the  first  and  second  kind  of  modulus  k,  and  F(4>.k')  and  E(4>,k')  ar<' 
the  incomplete  elliptic  integrals  of  the  second  kind  of  amplitude  and 
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modulus  (complementary)  k'. 

The  next  quantities  that  we  wish  to  determine  are  the  submergence, 
H-h,  and  the  spray  sheet  thickness,  6.  With  reference  again  to  Fig.  4, 
one  may  derive  the  submergence  by  combining  the  continuity  equation  and 
geometry  of  the  flow,  i.  e.  , 

H-h  =  d-e’,  (4.  17) 

where  e'  is  the  rise  in  free  surface  from  upstream  infinity  to  the  point 
N  at  the  distance  d  above  the  leading  edge  of  the  plate.  The  same  re¬ 
lations  give  the  spray  sheet  thickness. 

6  =  d  -  (e*  -e-Cha)  ,  (4  18) 

where  e  is  the  drop  along  the  lower  cavity  wall  from  the  trailing  edge  to 
downstream  infinity.  Comparing  the  above  expression  for  6  to  our  de¬ 
finition  of  g  given  in  Eq.  (3.4),  one  may  immediately  write 

g  =i(e' -e)  -  Ch  .  (4.19) 

We  shall  apply  our  present  theory  to  determine  e  and  e*. 

In  linearized  theory,  it  is  known  that  the  vertical  difference  be¬ 
tween  two  points  on  the  same  streamline  can  be  calculated  by  considering 
an  integral,  in  the  physical  plane,  of  the  form 

jvdx' 

where  values  of  v  are  taken  along  horizontal  lines  which  approximate  the 
actual  streamline.  In  our  problem,  all  free  surfaces  in  the  physical  plane 
are  mapped  onto  the  real  axis  of  the  £  -plane.  The  above  integral  can  then 
be  replaced  by  an  equivalent  one, 
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carried  along  the  real  4-axis.  This  enables  us  to  write  with  reference  to 


Fig.  3, 

00 

r% 

e' =  3^(0  “4  . 

(4.20) 

/»  b 

and 

(4.21) 

The  lower  limit,  n,  in  Eq.  (4.20)  is  the  image  of  N  in  the 
z-plane  and  has  to  be  determined  by  solving  a  transcendental  equation. 
Finally,  with  the  aid  of  the  mapping  function  in  Eq.  (4.  1  )  one  may  write 

e'^M,  (4.22) 


where  M  is  the  double  integral 
00 

( 

M  = 


f b+4 +  e-c ]  \/S+1f +b) \/(i+t)(b+t) d|  •  H-23) 


and 


(4.24) 


where  I  is  the  double  integral 


■f.h*  -rr]  \lEWEi 


dt  /  -t 

4-t  VTFTtyibTt) 


d  4 


(4.  25) 


All  that  now  remains  is  to  determine  n  and  to  evaluate  the  double  integrals 
M  and  I. 

i )  The  number  n  : 

For  the  calculation  of  n,  one  has  to  utilize  the  mapping  function  in 

Eq.  (4.  1).  Since  N  is  purely  imaginary  and  n  real,  the  substitution  of 

N  for  z  and  n  for  4  into  Eq.  (4.1)  gives 

.  ,  n  +  b  ,  n  -  c  „ 

b  log  — —  +  c  log  — —  0  , 
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or  (n/b  +  1  (n/c  -  1  )C  -  1  (4.26) 

Equation  (4.  26)  is  a  transcendental  one,  such  that  for  given  values  of  b 
and  c,  n  must  be  determined  by  iteration.  One  knows,  however,  that 
n  lies  in  the  interval  bounded  by  c  and  2  c. 
ii)  The  double  integral  M: 

The  inner  integral  in  Eq.  (4.  23)  can  also  be  expressed  in  terms  of 
complete  elliptic  integrals  of  the  firs*  and  third  kind.  If  we  further  re¬ 
place  the  integration  variable  £  by  1/s,  M  can  now  be  written  as  a 


single  integral, 
1 


1 


M  =  2K 


Cn  f-U  -r^-1  l-  d.  -  .  C"  p-  ♦ 

*o  i_s+k2  ~CS.  V  s(s+kl)  o  [s+k2 


A  (*,  k) 

o 


ds  , 

(4.27) 


where 


+  k2 


The  first  integral  in  Eq.  (4.  27)  can  easily  be  reduced  to  incomplete 
elliptic  integrals  and  M  now  can  be  expressed  as 


E {<p  ,  k' ) 


M  =  4K 


ck2  ( 1  +c) 


ck 


iZ 


JI(?  ,a  2  ,  k' )  + 
ck2  +  1  11  1  +  ck2 


F(?  ,k') 


_Tcnr_L 

•ils+k2  *'ci 


1  A(*,k) 


ds 


(4.  28! 


where  <p  =  sin  *  \  / — - —  <  ,  and  a  2  -  1  +  ck2>  1  • 

1  VUnk2  7  » 

The  remaining  integral  will  be  treated  numerically.  For  this  pur¬ 
pose  one  must  examine  the  behavior  of  A  /s  for  small  s.  By  definition, 


A(tt,k)  -  =-[  (E-K)F(*,k’)  +  KE(*,k')I  . 
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For  any  non-zero  k,  ^  approaches  zero  as  s  approaches  zero.  One 
may  therefore  apply  binomial  expansions  to  the  integrands  of  F(^,  k') 
and  k’)  and  integrate  term  by  term.  By  so  doing  the  asymptotic 

behavior  of  A^'I',  k)  as  ^  approaches  zero  is  found  to  be 


Ao(*,k)~|  E  *+  O(^)  .  (4.29) 

A  second  expansion  of  ^  is  then  carried  out  to  yield,  for  small  values 
of  s, 


* 


(4.  30) 


Hence,  by  keeping  only  the  first  term, 


A(^,  k)  ^  z  £ 

s  n  ^/s(s+k2 i 


(4.  5  1  ) 


We  have  therefore  found  that  the  term  2W  s  is  singular  at  s  =  0,  how  - 
ever,  the  singularity  is  of  the  order  of  one  half,  and  is  thus  integrable. 
In  order  to  facilitate  numerical  computation,  we  regroup  the  integrand 
in  the  following  form  ; 


1 


in 

1  c 

J 

O 

1  2  1  -  cs 

L  s  +  k  J 

A  (¥,  k) 


s 


ds 


■I 


n 


sfk^ 


-  cs 


2E 


AU.OI 


-  2E 


f 


h  J  L  Js(s4k2  ) 

ds 


ds 


o  L  s  +  k 


1  -  c  s 


\l Ms  +  k2  ) 


(4.  12) 


The  singular  behavior  of  the  integrand  is  thereby  subtracted  out  .  so  that 
the  first  integral  presents  no  difficulty  in  computation  and  the  second 
integral  can  be  integrated  in  closed  form. 
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Summarizing  the  results  for  M-integral,  we  have 


iii)  The  double  integral  I: 

The  treatment  of  the  integral  I  is  quite  similar  to  that  of  M. 

First,  the  inner  integral  is  expressed  in  terms  of  complete  elliptic  in¬ 
tegrals  to  reduce  I  to  the  single  integral, 

fb 

■  =  *4  [bTi  -rrrJ  r,-4«,-k»]  d*  ■  l4  54) 

where  P  =  sin 

The  whole  I  integral  will  be  treated  numerically.  As  in  treating  the  M 
integral,  we  shall  examine  first  the  behavior  of  the  quantity  (1-  A)/(b-  4) 
near  4  =  b.  With  the  aid  of  Legendre's  relation,  i.  e.  , 

EK'  +  E'k  -  KK’  =  j  , 

the  term,  1  -  A{{J,k)  vanishes  as  (3  approaches  j  ,  or  when  4  ap 
proaches  b.  One  can  show  further  that  in  the  neighborhood  of  b,  1  - 


fi-T 


0<  |5<£ 
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behaves  like  the  square  root  of  the  quantity  b  -  or  more  precisely, 


l-A  (P,k)~^  [|  -  Kk']  v/TTf 


(4.  35) 


We  therefore  apply  the  same  technique  as  we  did  in  the  treatment  of  \1 


integral  to  yield. 


where 


I  =  VI2 


(4.  36) 


b 

C  b 


V  V\  b^T  [l  -Ao0»k)]  ^ 

l 

=  [l-A(p.k)]-  -  Kk’] 

1 

+  4k'  [|  -  Kk’]  , 


(4.  36a) 


!2  *  -*e]  TTT  [‘  -Ao<P’k>]  ■ 


(4.  36b) 


The  numerical  computations  are  described  in  detail  in  the  appendix. 
All  numerical  integrations  are  carried  out  on  the  IBM  7090  computer  at 
the  computing  center  at  the  California  Institute  of  Technology. 

5.  Limiting  cases 

The  present  solution  contains  several  geometrical  parameters  such 
as  h,l  ,  etc.  When  these  parameters  approach  zero  or  infinity  several 
interesting  but  simpler  flows  are  obtained  Recovery  of  results  for  simple 
flows  from  solutions  of  a  more  general  one  often  serves  as  an  indication 
of  the  correctness  of  the  general  solution  and  also  provides  one  certain  u>» 
ful  information  about  the  behavior  of  the  solution  near  such  limits  We 
shall  consider  several  cases  The  parameters  to  be  varied  are  the  dis- 


la 


tances  h  and  d,  and  the  ratios  C,  /d  and  C,  /h-  In  all  cases,  we  are 

n  h 

interested  in  the  behavior  of  the  slope  of  the  lift  coefficient  C.  The 

Hi 

length  of  the  wall  I  is  equal  to  positive  infinity,  i.  e.  a  completely  rigid 
bottom  for  the  first  four  cases  and  negative  infinity  (or  a  free  jet)  for  the 
last  one. 

All  cases  are  for  vanishing  angle  of  attack  so  that  we  may  put 

6  =  d. 

Case  1:  The  infinitely  deep  stream  (h  -*oo)  . 

Given  a  plate  of  finite  length  C^<  the  limit  Cy/h~*  0  implies  that 
h  or  b  approach  infinity,  or  that  the  stream  becomes  infinitely  deep. 

C.  as  a  function  of  C,  /d  is  then  given  by  the  simultaneous  equations 

a  h 

Ch/d  =  J  "  -f  loS  (1  +  T/d) 


c^a- 


as  may  be  found  from  Eqs.  (4.  9)  and  (4.  15).  These  results  agree  with 

(2) 

those  of  Schot  Again  it  should  be  noticed  that  as  the  submergence 
approaches  zero,  tends  to  the  limit  it,  and  for  the  infinite  fluid 

case,  approaches  w/2  . 

Case  2:  Planing  on  a  stream  of  finite  depth  (  d  0  ). 

Since  d,  the  spray  sheet  thickness,  vanishes,  this  case  corres¬ 
ponds  to  that  of  planing.  We  obtain  for  this  limit 


4 

C,  IT 
n 


K2 
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near  d  =  0.  Introducing  a  new  parameter  k^  ,  sui  Is  that 

!  -  k' 


k  -  _ 

1  1  +  k’ 


Reference  (7)  gives  the  following  expansion  ol  K 

k  =-j(i  +  k4)  ,  \ ;  *' S2} 


where  F  is  the  usual  hypergeometric  function. 

Now  we  let  the  depth  of  the  stream  become  large  and  the  leading 
terms  of  the  following  expansions  are  readily  obtained  : 


and 


k*  =  4‘b  (‘-i)  *°(b3)' 

=  +2b  ♦  3?  ^ 

W]- 


IT 

L  =  cT 

a  h 


,  1  11  1  A 

1  +  -rr-  +  -r-=r  -  +  0 

2b  32  ^2 


From  the  mapping  relation,  we  have 

tt  C, 


*('  *^)  •  °(?) 


and  if  the  above  relations  are  combined  for  small  ratios  of  C^/li  ,  \w 


have  the  approximate  result 


CT  -ir  1  i  IM  (,2H0h  (C. 

hi  L  h 


The  first  term  is  the  result  » »f  planing  with  inlmite  depth  at  inlinite 

(5) 


Fronde  number.  y'f  The  total  expression  t < > r  C  ignoring  terms  ler>s 

a 

than  l)*2  is  very  close  to  the  result  <>t  J).i v.  son  ( 1  )  winch  gives 

cL  -  [i  •  O.IOJ  fCj./i.e] 
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Case  3:  Flow  past  a  cavitating  plate  near  a  solid  wall  (d-»oe). 

This  limit  corresponds  to  the  case  of  infinite  submergence.  The 

plate,  however,  is  held  at  an  arbitrary  distance  from  the  wall.  C .  as 

a 

a  function  of  the  distance  h  (or  b)  is  now  governed  by  the  simultaneous 
equations 


and 


4b 

h 


(K  -  E)2 


Ch=  -E  1°S  >• 

We  will  now  obtain  the  first  order  wall  effect  by  letting  h  become 
large.  We  obtain  after  a  little  manipulation  the  approximate  result 

CL  =  J  f  1  +  °-208  \l  Ch/h 

a 

The  constant  term  again  is  a  well  known  result  of  infinite  fluid  flow. 
Case  4:  The  case  of  very  large  chord  (C^  "•’oo  ). 

There  are  two  simple  limiting  situations:  One  of  zero  submergence 
and  one  of  infinite  submergence.  We  will  consider  them  separately. 

i)  Zero  submergence.  We  may  set  c  in  the  mapping  function 
to  be  zero  to  obtain 


The  slope  of  lift  coefficient  now  becomes 


exp 


ttC, 


C 


L_ 


4 

V 


K* 


for  this  case.  Now  we  will  let  the  chord  become  large  compared  to  the 
depth,  so  that  C^/h  approaches  infinity  and  k'  approaches  zero.  The 
K  function  has  then  the  well  known  behavior 

K  log  4/k\ 


1 1 


As  b  —  1 , 


We  get  the  same  result  as  the  previous  case. 

Case  5:  The  free  let  (f  -x)  . 

The  entire  flow  in  this  cast*  be  onies  a  free  jet  The  parameter 

o'  in  Eq.  (4.  10)  approaches  plus  infinity  as  the  bottom  wall  v  ir.mh, 

and  the  asymptotic  angles  of  inclination.  0^  and  0  given  m  Eqs  t4  7  ) 

and  (4.8),  can  be  integrated  in  closed  form  to  gi\e  the  follow  my  result 

for  C  ^  with  the  use  of  Eq.  (  T  5) 

a 
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This  result  and  that  of  mapping  function  (Eq.  4.  9)  give  the  com¬ 
plete  solution  for  the  cavitating  flow  past  a  flat  plate  hydrofoil  in  a  free 
jet.  As  a  further  special  case  to  illustrate  the  effect  of  a  lower  free  jet 
let  us  take  the  submergence  to  be  infinite.  Then,  c  approaches  infinity 
and  the  mapping  function  Eq.  (4.  9)  becomes 

ch/h=- iiog  (l-p)-  ' 

When  the  hydrofoil  approaches  the  bottom  closely,  b  approaches  unity 
and  the  following  approximate  formula  for  CT  can  be  quickly  establish- 

CV 

valid  for  large  values  of  C^/h.  Thus,  the  lift  slope  vanishes  when  the 
hydrofoil  is  on  the  lower  free  surface  and  grows  linearly  with  depth 
for  small  submergences. 

6.  Numerical  Results 

The  majority  of  the  numerical  results  presented  herein  are  for 
the  case  of  the  infinitely  long  bottom  wall.  To  bring  out  the  salient 
features  of  this  more  common  situation,  we  have  chosen  several  values 
of  C^/h,  within  the  usual  range  of  experimental  interest,  and  calculated 
the  lift  slope,  C.  ,  the  spray  sheet  contraction,  6/d,  as  functions  of 
the  chord  spray-sheet  thickness  ratio,  C^/d,  and  the  dimensionless  sub¬ 
mergence  ratio,  (H-h)/C^.  These  are  shown  in  Figs.  5,  6  and  7  res¬ 
pectively. 
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An  inspection  of  Fig.  5  brings  out  the  fact  that  the  effect  of  the 
bottom  and  that  of  the  free  surface  are  in  the  same  direction;  namely, 
to  increase  the  lift  slope.  Moreover,  the  influence  of  both  the  bottom 
and  the  free  surface  extends  to  many  chords  and  decreases  even  then 

at  a  slow  rate.  For  example,  the  fully  submerged  limit  (C^  =  tr/2)  is 

a 

not  reached  even  with  a  submergence  of  100  chords;  and  the  presence  of 
the  bottom  is  still  in  evidence  with  ten  chords  of  depth  beneath  the  foil. 

Interestingly,  the  spray  sheet  contraction  6/d  show's  only  a 
minor  dependence  on  channel  depth  ratio,  C^/h,  but  it  is  affected  to  a 
marked  degree  by  the  angle  of  attack.  This  is,  of  course,  just  the  effect 
of  g  as  indicated  in  Eq.  (3.  4).  In  fact,  it  would  have  been  perhaps  mor< 
illustrative  to  plot  the  ratio  of  g/d  rather  than  that  of  6/d  as  is  done 
in  Fig.  b.  In  any  event,  even  at  the  moderate  angle  of  six  degrees,  the 
spray  sheet  can  become  significantly  less  than  distance  d.  We  may  add 
that  previous  linearized  theories  have  not  distinguished  between  these 
qua  ntities. 

The  submergence  of  the  hydrofoil  beneath  the  undisturbed  sur- 
face  is  given  in  Fig  7  for  various  ratios  of  C^/d  and  chord-depth 
ratio  C^/h.  The  submergence,  like  the  spray  sheet  will  depend  upon 
the  angle  of  attack.  At  vanishingly  sum  11  angles,  the  subnu  rgence  be¬ 
comes  equal  to  the  quantity  d.  but  even  at  small  angles  it  is  seen  to 
become  appreciably  less  than  this  The  submergence  ratios  in  Fig.  7 
are  calculated  for  a  -  6°.  The  limiting  case  \  0  is  also  shown.  It 

is  not  surprising  that  for  small  values  of  C^/li  and  large  values  of 
Cj  /d  that  the  submergence  tan  become  negative  even  leading  to  the 
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situation  in  which  the  entire  hydrofoil  lies  above  the  undisturbed  water 
level.  This  result  is  well-known  from  the  previously  mentioned  work  of 
Green. 

For  the  finite-length  bottom  calculations,  we  chose  the  case 
where  the  plate  is  in  midstream,  that  is  b  =  c,  and  C^/h=  0.4,  be¬ 
cause  this  situation  corresponds  to  the  usual  experimental  conditions 
in  the  Free  Surface  Water  Tunnel  at  C.I.  T.  Figure  8  shows  C 
versus  I  / for  this  case.  When  the  bottom  extends  to  approximately 
six  chords  downstream  from  the  leading  edge  of  the  flat  plate  the  free 
jet  effect  becomes  negligible,  so  that  the  lift  slope  is  essentially  the 
same  as  that  calculated  for  the  infinite  bottom  case.  Similarly,  at 
approximately  six  chords  upstream  the  lift  slope  approaches  that  of  the 
completely  free  jet.  This  result  shows  that  the  length  of  bottom  required 
to  establish  the  infinite  bottom  effect  is  in  reality  only  a  few  chords  -  say 
about  three  downstream  of  the  hydrofoil.  Interestingly,  nearly  the  in¬ 
finite  fluid  case  is  achieved  when  the  bottom  stops  just  underneath  the 
leading  edge  of  the  hydrofoil.  Seemingly,  the  free-surface  effect  and 
the  bottom  effect  are  then  just  about  cancelled  out. 

The  details  of  the  present  calculations  are  discussed  in  the 
Appendix  and  the  intermediate  parameters  needed  for  the  calculations 
and  for  calculations  at  other  angles  of  attack  than  given  so  far  are  tabu¬ 
lated  in  Tables  1  and  2.  The  geometrical  parameters  b  and  c  and  the 
dependence  of  J  on  the  quantities  are  graphed  in  Fig.  A-l  of  the  Appen¬ 


dix. 
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7.  Discussion 

It  would  be  desirable  to  be  able  to  present  experimental  data  con¬ 
firming  the  theory  just  outlined.  The  ubiquitous  presence  of  gravity, 
however,  makes  all  such  hopes  faint;  zero  cavitation  number  - -resulting 
in  physically  unrealistic  infinite  jets- -is  not  possible  in  a  gravitational 
environment.  The  indicated  lift-slope  values  of  Fig.  5  for  geometries 
similar  to  those  used  by  Dawson  and  Bate  are  far  higher  than  those 
extrapolated  from  their  experiments  which  were  conducted  at  non-zero 
cavitation  numbers.  Undoubtedly  the  effects  of  the  free  surface  and  the 
rigid  bottom  are  greatly  augmented  for  zero  cavitation  number,  in  fully 

wetted  flow,  a  relatively  slight  effect  on  for  the  geometry  of  Fig.  8 

a 

would  be  expected,  for  example.  Nearer  the  free  surface  it  becomes 
possible  to  approach  more  closely  the  conditions  of  zero  cavitation  number 
so  that  the  conditions  of  the  theory  will  be  more  accurately  met  there. 

To  investigate  this  point  a  flat  plate  hydrofoil  was  mounted  in  a 
two-dimensional,  parallel  wall  insert  six  inches  ide  in  the  Free  Surface 
Water  Tunnel  similar  in  type  to  the  one  used  by  Dawson  and  Bate  Two 
photographs  of  the  hydrofoil  and  resulting  spray  sheet  are  shown  in  Fig  9 
In  these  photographs  three  water  surfaces  can  be  seen  the  level  exterior  to 
the  two-dimensional  insert,  the  lower  and  upper  boundaries  of  the  cavity  . 
and  finally  the  top  of  the  spray  sheet  itself.  The  curving,  exterior  water 
surface  just  above  the  hydrofoil  should  be  disregarded.  It  should  be  men¬ 
tioned  that  the  effects  of  friction  on  the  inside*  of  the  channel  result  in  a 
very  gradual  curvilinear  flow,  the  surface  of  which  is  slightly  uu  lined  to 


the  horizontal.  Angles  of  attack  were  determined  relative  to  tins  undisturb 
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ed  flov.\  It  was  not  possible,  however,  to  make  sufficiently  precise  mea¬ 
surements  of  the  effective  angles  of  attack  to  permit  C ^  to  be  calculated. 
Measurements  of  the  distance  d  and  submergence  below  the  undisturbed 
surface  were,  however,  made  and  these  are  given  in  Fig.  10  for  several 
angles  of  attack.  Also  shown  are  the  theoretical  values  for  the  indicated 
angles,  calculated  for  the  geometry  of  the  experiment.  Data  were  taken 
at  two  speeds,  12  ft. /sec.  and  18  ft. /sec.  The  chord  of  the  hydrofoil  was 
four  inches. 

It  is  apparent  from  Fig.  10  that  the  observed  submergence  is  al¬ 
ways  greater  than  that  calculated,  and  that  the  discrepancy  increases  with 
angle  of  attack.  There  is  no  perceptible  difference  in  behavior  with  ve¬ 
locity  (or  Froude  number)  although  one  may  observe  in  Fig.  9  that  there 
is  a  substantial  gravity  effect  on  the  spray  sheet  causing  it  to  follow  a 
curved  trajectory.  The  actual  magnitude  of  the  difference  in  the  observed 
and  calculated  value  at  6°  angle  of  attack  amounts  to  about  0.  1  chord  or  in 
the  case  of  the  present  experiment  nearly  0.  4  inches.  This  is  too  large  for 
experimental  error.  In  these  spray  sheet  measurements  the  cavity  was 
vented  to  atmosphere.  Nevertheless,  there  was  still  a  hydrostatic  gradient. 

We  are  at  somewhat  of  a  loss  for  a  quantitative  explanation  of  the 
lack  of  agreement  between  the  theory  and  experiment  insofar  as  :he  spray 
sheet  measurements  go.  The  experiments  do,  however,  follow  the  pre¬ 
dicted  trends,  even  exhibiting  the  negative  submergence  although  to  a 
smaller  magnitude.  Gravity  would  have  such  an  effect  as  it  tends  to  pre¬ 
vent  the  free  surface  from  rising.  At  the  same  time,  the  observations 
suffer  because  the  measurements  were  made  in  a  rather  narrow  channel 


11 


(six  inches).  The  sidewall  boundary  layers  in  the  channel  at  the  site  of 

(9) 

the  hydrofoil  amounted  to  about  one -half  inch  in  thickness  and  there 
was  also  a  small  but  undetermined  curvature  of  the  undisturbed  free 
surface.  The  slope  of  this  surface  at  the  point  of  hydrofoil  immersion 
was  three  degrees.  With  this  combination  of  circumstances,  one  should, 
perhaps,  not  expect  better  agreement  with  a  linear  theory  at  zero  cavita¬ 
tion  number. 

In  closing,  it  is  worthwhile  to  mention  that  it  is  easy  to  extend 
the  present  calculations  to  account  for  arbitrary  hydrofoil  shape.  Further 
more,  it  is  not  much  harder  to  carry  the  problem  through  for  non- zero 
cavitation  numbers  provided  one  of  the  simpler  "wake"  models  is  adopted 
as  the  prescription  for  the  flow. 
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APPENDIX.  CALCULATION  PROCEDURES 

A.  Infinitely  Long  Bottom  Wall 

The  first  step  in  the  calculation  of  the  desired  quantities  C-t 

L- 

a 

(H-h)/C,  ,  6/d,  and  C ,/d  is  to  determine  values  of  b  and  c  for  each 

n  h 

specified  C^/h.  This  entails  solving  by  iteration  the  transcendental 
equation  given  in  Eq.  (4.  9).  We  should  point  out  that  once  a  value  of 
C^/h  is  prescribed,  the  values  of  b  lie  within  a  particular  range  and 
no  solution  outside  that  range  is  possible.  The  solid  lines  in  Fig.  A-l 
are  curves  of  constant  C^/h  in  the  b-c  plane.  Notice  that  as  C^/h 
approaches  unity,  the  range  of  b  for  which  a  solution  exists  becomes 
very  narrow.  These  curves  may  be  used,  along  with  additional  values 
given  in  Table  1,  to  obtain  first  approximations  of  b  and  c  for  values 
of  C^/h  not  tabulated. 

Once  a  set  of  b's  and  c’s  are  established,  for  the  infinitely 

long  bottom  wall  case,  it  is  now  easy  to  solve  Eq.  (4.  16a)  or  (4.  16b)  for 

J  and  then  obtain  from  Eq.  (4.  15).  If  J  is  calculated  on  a  high 

a 

speed  computer,  Eq.  (4.  16a)  is  a  convenient  form  to  use  as  a  computer 
method  for  evaluating  Il(a2  ,  k)  is  given  in  Ref.  (8).  Or,  the  other  hand 
Eq.  (4.  16b)  is  a  more  suitable  form  to  use  when  J  is  calculated  on  a 

desk  calculator  since  K  and  A(4».k)  are  tabulated  in  Ref  (7)  Tin- 

o 

dashed  lines  in  Fig.  A-l  are  lines  of  constant  J  in  the  b-c  plane. 

The  submergence  ratio  (H-h)/C^  is  calculated  from  Eqs  (4  17). 
(4.22)  i4  26b)  and  (4.  33).  I’he  exact  part  of  the  M  integral  us  given 
in  Eq.  (4  3})  is  easily  evaluated  by  the  methods  described  in  Ref  (8t 
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Since  the  integral  part  of  M  is  always  less  than  two  percent  of  the  total 
M,  (at  least  for  the  range  of  cases  that  we  calculated)  the  numerical 
integration  need  not  be  done  to  a  great  degree  of  accuracy.  When  the 
integral  is  divided  into  three  appropriately  chosen  intervals  it  converges 
very  rapidly  to  the  desired  accuracy  of  a  total  of  five  significant  figures 
in  M  or  at  most  three  figures  in  the  integral  part. 

The  spray-sheet  contraction,  6/d  requires  the  calculation  of 
the  I  integral  as  well  as  the  M  integral.  The  I  integral  as  given  in 
Eq.  (4.  25)  is  divided  into  two  integrals  1^  and  I  ,  where  1^,  the 
more  lengthy  calculation,  is  independent  of  c  and  hence  independent  of 

ch/h- 

B.  Finite  Bottom 


In  order  to  calcu’ate  the  lift  slope  for  this  case  the  parameter  b' 
must  be  evaluated  from  Eq.  (4.  10)  for  each  prescribed  t  (or  B")  and 
any  set  of  C^/h,  b,  and  c  calculated  as  outlined  above.  We  chose  to 
do  only  one  case  where  b  =  c  and  hence  Eq.  (4.  10)  could  be  solved 
directly  for  b'  rather  than  by  iteration 

The  lift  slope,  C.  is  now  readily  obtained  from  Eq  (3.  5)  in 

a 

which 


and 


where 


0 

i 


a  \J( c+  1  )(b'  +  c  ) 


a  .  a,  /(b-D(b'-b) 

®  -  f  v — b — 


f1  r  t  dt 

J .  V  (T  - 1  )(b-  ^7)  t  ♦  c 


\/<)  -t)(b'-u  b-t 


J  and  j’  can  be  reduced  to  elliptic  integrals  with  the  following  results 


-k  ["  ("f-  •<  )-K(k  .] 


or  in  alternate  form 


“yil-cTIbSc)  V*'-  k  >-t? 


Zj£ 

b’  +  c 


and  for  J*  we  find 


j-  2k  jn(-i  ,  k  )  -  K(k  )  j 


J'  -  k  k 


1  -  A  { 0 ,  k  ) 

V(k*  -  k*  2 )( 1  -  k2  ) 


where 


o  s  m 
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‘V  i  -  k"2 
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Figure  1  .  Sketch  showing  a  cavitating  flat  plate  in  a  stream  of  finite 
depth  with  the  lower  surface  partly  free  and  partly  solid. 
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Diagram  of  the  physical  z-plane  showing  boundary 
conditions  for  the  linearized  flow. 
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Figure  4.  Detinitiun  sketch  of  a  cavilating  flat  plate  in  a  stream  of  finite  depth. 
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Figure  5.  The  lift  slope  versus  C^/d  for  various  values  C^/ 

a 

for  a  flat  plate  hydrofoil  in  a  stream  of  finite  depth. 


Figure  7.  Submergence  of  the  plate  versus  C^/d  l"r  several 
values  of  C^/h  at  an  angle  ol  attack  of  6* 
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Figure  8.  The  lift  slope  CT  versus  various  length  of  bottom  wall  with  the  hydrofoil 

midway  between  top  and  bottom  surfaces  and  for  the  chord-depth  ratio  of  0.4. 


Figure  9.  Two  photographs  of  the  planing  flat  plate  at  an  angle  of  attack 
of  about  twenty  degrees  at  a  velocity  of  about  12  ft.  per  see. 


EXPERIMENTAL  DATA 


Submergence  of  the  plate  versus  spray  sheet  thickness  at 
the  leading  edge  for  various  angles  of  attack.  The  solid 
lines  are  the  theory  evaluated  for  a  value  o  f  C,  /h  =  0.  2. 
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9.2912  X  10  * 

59  41) 

1  040 

1-10041 

1.0004 

1  02)5 

1 .33590  X  10 

2  9124  X  10  * 

1.045 

7.52005  X  10’* 

2.3120 

1  024 

6  99328 

4  6095  X  10'* 

29  596 

l  070 

3. 371 14  X  10'* 

2.4545 

1  0245 

6.29088 

6  3261  X  10  * 

26  216 

1  074 

3.04700  X  Ur* 

2.04.44 

1  025 

4  79589 

7  9672  X  10 

24  no 

1.000 

2  04231  X  10'* 

3  1005 

1  0)0 

1.21719 

2  1091 

17  770 

1.004 

2.20200  X  10'* 

3.4210 

1  0)5 

5.73621  X  1C'* 

3  0112 

16.190 

1.090 

1  44145  X  10  '* 

34171 

14  467 

4  277 

1  040 

3. 11688  X  10'* 

3.6653 

15  414 

MOO 

0.41313  X  10** 

1  0527 

12.700 

4.204 

1  045 

1.73797  X  10'* 

4.1613 

14  929 

1.105 

5.01004  X  10-* 

4.1005 

1  050 

9  22167  X  10-* 

4.5542 

14.540 

1.110 

3.74431  X  10'* 

4  2406 

1  055 

4  1  3876  X  10'* 

4  0433 

14  295 

1.115 

2.07742  X  10'* 

4.3785 

1  0575 

2.17190X  1C'» 

5  0  36) 

14  161 

1  1175 

1.1*704  X  10'* 

4.4438 

1  060 

1.02074  X  10  * 

5  1935 

14  019 

1.120 

7.00035  X  lO* 

4.5242 

12.274 

4  294 

1  0605 

7  98235  X  10'* 

*  2273 

13  944 

MSI 

5.00404  X  10'* 

4  5305 

1  0610 

5.92171  X  10'* 

5  2630 

u  954 

1.122 

4 . 000  3  3  X  lO* 

4  5957 

1  0615 

4.03433  X  10'* 

5.  *315 

i  3  914 

M2) 

2  30202  X  10'* 

4.4348 

l  062 

2  33771  X  10'* 

5  3*47 

13.877 

1  0625 

8.71556  X  10  '* 

5  3986 

1  3  424 

TABLE  2 

t»cl.  182302  Ch/h  0  4 


b' 

J# 

J' 

CL 

a 

L/C, 

1  143 

1.3606 

9  2919 

2  05 

5  )6 

1  145 

1  3562 

9  2542 

2  04 

4  29 

1  190 

1 . 3452 

9  1620 

2  0) 

)  45 

1  200 

1  3243 

4  9663 

2  01 

2  79 

1  220 

1  2156 

4  6656 

197 

2  14 

1  250 

1  2  344 

3  2471 

1  91 

\  70 

i 

1  1627 

7  67  34 

1  44 

l  25 

1  400 

J  0533 

6  4260 

1  74 

0  72 

1  500 

0,9714 

6  2166 

1  66 

C  )9 

2  00 

0  7417 

4  5469 

1  49 

0  50 

3  00 

0  5509 

3  3292 

1  )8 

-1.35 

4  00 

0  4544 

2  7450 

1 

1  87 

6  00 

0  3609 

2  1440 

1  30 

-2  55 

4  00 

0  3073 

1  4174 

1  29 

■  )  0) 

10.00 

0  2722 

l  6074 

1.28 

<  19 

15  00 

0. 2  1  ?4 

1  29.12 

1  27 

•4  04 

20  00 

G  1644 

1  1106 

1  26 

4  50 

30  00 

0  15)2 

0  9001 

1  254 

-5  >5 

50.00 

0  1141 

0  69)1 

1  250 

5  96 

100  00 

0  04  32 

0  4440 

1  25 

-7  06 

